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ABSTRACT 

We present multi-epoch, total intensity, high-resolution images of 43GHz, v=l, 
J=l-0 SiO maser emission toward the Mira variable R Cas. In total we have 23 epochs 
of data for R Cas at approximate monthly intervals over an optical pulsation phase 
range of <fi = 0.158 to <fi = 1.78. These maps show a ring-like distribution of the maser 
features in a shell, which is assumed to be centred on the star at a radius of 1.6 — > 2.3 
times the radius of star, i?*. It is clear from these images that the maser emission is 
significantly extended around the star. At some epochs a faint outer arc can be seen 
at 2.2 R+. The intensity of the emission waxes and wanes during the stellar phase. 
Some maser features are seen infalling as well as outflowi ng. We have made initial 
comparisons of our data with models bv lGrav et. al.l (|2009t ). 

Key words: maser - techniques: interferometric - star: AGB stars and post- AGB - 
star: individual: R Cas . 



1 INTRODUCTION 

The Asymptotic Giant Branch (AGB) is the last evolution- 
ary stage of low- and intermediate-mass stars that are nu- 
clear powered. This phase of evolution is characterised by 
nuclear burning of hydrogen and helium in thin shells above 
the electron degenerate core of carbon and oxygen. These 
stars are very large (R~l AU), very cool (2000-3000 K) 
and have abundant, mol ecular, dusty winds (mass loss rate 
10" 7 - 10" 5 M yr _1 jHerwidl2005l ). 

Any AGB star in the end will lose mass in the form of 
a slow wind at a rate that will significantly affect the mass 
of the star. This produces a circumstellar envelope of es- 
caping dust and gas particles. When the envelope becomes 
separated from the star, the mass-loss sometimes ceases in a 
short period of time and the temperature drops as the wind 
flows away from the star, from about 3130 K at the pho- 
tosphere to « 740 K at 5.5R* (|Grav et. al.l 120091 ') ; cooling 
takes place by line radiation from various molecules, espe- 
cially H2O. Some molecular species are formed in an equi- 
librium process deep in the atmosphere and are destroyed in 
the outer parts of the outflow by interstellar ultraviolet ra- 
diation (H2, CO, H2O). Other species, for example SiO, are 
depleted due to conden sation on dust particles at a few stel- 
lar radii (|Habinell 199(1 . The observations show that almost 
all AGB stars are pulsating variables, the pulsations pro- 
duce waves which steepen up into shock waves in the large 
density gradient of the stellar atmosphere. These pulsational 
shocks are well-know n to influence the dust formation zone 
|Frevtag et. al.l [20081 ') . Stellar pulsations cause waves which 



travel into the thinner outer layers of the atmosphere and 
develop into shocks with temperature and density variations 
strong enough to trigger intermittent dust formation. Large 
opacities allow the dust grains to capture the stellar radi- 
ation efficiently, leading to an outward acceleration of the 
dust particles which is transmitted to the gas by collisions, 
causing a slow but dense outflow from the star. The mass 
loss controls the termination of the AGB phase. When the 
envelo pe becomes optically t hin, the star leaves the AGB 
phase (jLattanzio et. allll998l ). Shock waves can produce the 
necessary physical condition (SiO abundance, kinetic tem- 
perature and density ) for collisional pumping of SiO maser 
|Elitzur et. al| [l980). Their positions is also consistent with 
that of optically thick molecular layers responsibl e for ra- 
diative pumping of these masers l|Grav et. al.ll2009l ). 

The gas in the circumstellar envelope is largely molec- 
ular with a composition that reflects the stellar chemistry. 
The molecules are located in different parts in the envelope 
in shells of different radii and widths, depending on the for- 
mation p rocess and thei r resistance to shocks and ultraviolet 
photons l|01ofssonlll999l ). 

The gas outflow of AGB star is normally considered 
to be spherically symmetric. Asymmetric structures are ob- 
served in AGB stars, for exa mple in the radio photo spheres 
of Mira, W Hya an d R Leo dReid fc Mentenll2007l ) and in 
optical photosphere (|Karovska 19971 ). 

R Cassiopeia is an oxygen-rich AGB star which is classi- 
fied as an M-type Mira- variable. Its optical brightness varies 
from magnitude +4.7 to +13.5 with a period of 430 days and 
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its mass is about 1.2 Mg. IVlemmings et all (|2003l ) used as- 
trometric VLBI to measure a distance of 176^4?, P c with a 
proper motion of (85.5 ±0.8, 17.5 ±0.7) mas yr" 1 in R. A. 
and Dec, respectively. Various estimates of the stellar veloc- 
ity V* appear in the literature; we tabulate these in Table [All 
and adopt the mean value, ±24 km s , standard deviation 
2 km s" 1 . 

Because the maser emission is so bright, it allows 
the inner circumstellar envelope to be resolved on sub- 
milliarcsecond scales. The SiO masers are found in a region 
close to the star, within the dust formation zone, located 
at a radius of about 2 to 4 stellar radius (R*). They are 
observed as clumpy, partial rings centred on the star; the 
masers fit a model of a hollow, expanding ellipsoid as ob- 
served in TX Cam and U Her. The masing action needs a 
long path through a relatively constant tangential velocity 
to produce the necessary gain to be visible, this condition is 
often met for a line tangent to the masing region with low 
turbulence because the turbulence would reduce the path 
length at a given velocity limitin g the gain and hence the 
brightness (|Diamond et. al. lll994T ). 

The H2O maser is located in the inner parts of the cir- 
cumstellar shell between w 10 to 20 R* where the tempera- 
ture r anges from 300 K to 1000 K and d ensities are 10 7 -10 9 
cm i Benson fc Little- MarenirJ[l996h . The OH masers are 
found i n the region furthest from the star as described by 
ICohenl (|l989l ). 



2 OBSERVATIONS AND DATA REDUCTION 

The SiO masers around R Cas were observed as part of a 
more extensive programme of VLBAQ monitoring of other 
stars. In total we have 23 epochs of data for R Cas. Data 
were recorded at each VLBA antenna in dual-circular po- 
larization in two 4 MHz windows, each digitally sampled at 
the full Nyquist rate of 8 Mbps in 1-bit quantization. The 
lower spectral window was centred at a fixed frequency cor- 
responding to the v=l, J=l-0 SiO transition, at an assumed 
rest frequency of 43.12207 GHz and a systemic velocity Vlsr 
= ±24 km s" 1 . R Cas was observed for three 45-minutes 
periods evenly spread over the 8 hour duration of the run. 
Adjacent to each R Cas observation, 5 minutes was spent 
observing the continuum calibrator 0359±509 at the same 
frequency as R Cas. The data were correlated at the VLBA 
correlator in Socorro, NM. The correlator accumulation in- 
terval was set to 2.88 seconds. All polarization correlation 
products (RR, RL, LR, LL) were formed. This configura- 
tion produced auto- and cross-power spectra in each 4 MHz 
baseband with a nominal velocity spacing of ~ 0.2 km.s -1 . 
Table [1] shows the epoch code, the date of observation and 
the corresponding optical phase. 

We reduced the data using the standard ap- 
proach to VLBI spectroscopy within the NRAO AIPS 
package (http://www.aips.nrao.edu/cook.html). We pro- 
cessed the visibility data using the the semi-automated 
spectral-line polarization calibration pipeline described 
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and further developed by 
iKemball fc Diamond! (|l997f i. We summarise the main steps 
here. We started by editing the data. We derived correc- 
tions for the phase-rate and the delay using the point-like 
calibration source 0359±509. The Kitt Peak (KP) antenna 
has good stability in both pointing and receiver gain so we 
used it as the reference antenna for all relevant stages at all 
epochs. We also scaled the KP autocorrelation spectra for 
0359±509 with respect to the system temperature in order 
to derive the amplitu de scale for all cross-power data (see 
IKemball et~al 1 , 200$ ) for details). We also used 0359±509 
to correct for the phase response across the bandpass. We 
applied the delay, rate, amplitude and bandpass calibration 
to R Cas. 

At each epoch, we selected a reference group of strong 
channels from R Cas and derived the residual rate correc- 
tions and then phase and amplitude corrections, which were 
applied to all channels. We took care to preserve the low 
Stokes V signature during amplitude calibration. We also 
corrected for parallactic angle rotation. Additional polar- 
ization calibration and results will be presented in a sepa- 
rate paper. After all calibration, we produced image cubes 
for R Cas at each epoch. We imaged 56 channels around 
the peak, each 0.217 km.s -1 wide, image size 102.4 x 102.4 
or 204.8 x 204.8 mas, in order to enclose a greater extent 
than any expected emission. The naturally-weighted restor- 
ing beam varied from epoch to epoch by a few mas and by a 
few degrees in position angle, but the close similarities show 
that the visibility plane coverage was consistent. The aver- 
age beam size was (0.039 ± 0.007) x (0.020 ± 0.008) mas 2 at 
average position angle (39 ± 17)°. The thermal noise <r rms in 
quiet channels was 30 mjy beam" 1 . 

We used the AIPS task SAD to fit 2-dimensional Gaus- 
sian components to each clump of emissi on brighter than 
5_£rma in each channel, using the method of lGonidakis et. all 
(2010). We grouped components if their position in adjacent 
channels agreed within a beamwidth ~ 1 mas, rejecting iso- 
lated components. We term these groups 'features' and we 
found the error-weighted mean position and flux-weighted 
mean velocity of each feature. Since there was no external 
phase-referencing, all positions were initially obtained rela- 
tive to the brightest component at each epoch. We aligned 
the emissi on at different epoc hs following the method de- 
scribed bv lCotton et. alj (|2004l L For each candidate position 
we found the flux-weighted mean radius of all SiO maser 
components and then the inner and outer radii which en- 
closed 90 percent of the total emission. The centre of ex- 
pansion is defined as the position giving the narrowest ring 
enclosing 90 percent of the emission. We then inspected the 
changes in component distribution relative to each centre. 
Some series of components formed very distinctive patterns 
which could be used to trace the emission from one epoch to 
the next. In most cases, alignment using the fitted centres 
of expansion gave very small and approximately symmetric 
shifts. However, at some epochs the maser distribution was 
very lop-sided, the shells were very poorly defined and we 
adjusted the apparent centres by eye to minimize the posi- 
tion changes of distinctive features. The estimated radii are 
given in Tabled We estimated the errors arising from: 

(a) Signal-to-noise ratio based fitting errors from SAD; 

(b) Difference between finding the centre by eye and the 
automated fitting. 
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Epoch Code 


Observing Date 


Optical Phase 


RD62 A 




0.158 


RD62R 


1 QQQ_Ort-1 5 


0.241 


RD62C 


1 QQQ-Nov-1 4 


0.310 




1 QQQ Tier 1 Q 

J. ffj JJCL _L £7 


0.390 


RD62F 


2000- Tan-1 5 

iUWU U till _1_ 


0.452 


RD62F 


2nnn_Fph-i 4 


0.521 




9000 A/Tar 1 7 


0.595 


RD62H 




0.675 


BD62I 


9000 Mflv 91 


0.744 


RDfiQA 

1 J 1 -> V J £7 . i 


2DnD- Tnn-22 


0.818 


RD69R 


2nflD- Tnl-1 fi 

v / v / V » • y i i i i v ' 


0.873 






0.956 


RD69D 


2000-Sep-21 


1.027 


BD69E 


2000-Oct-19 


1.091 


BD69F 


2000-Nov-17 


1.158 


BD69G 


2000-Dec-20 


1.234 


BD69H 


2001-Jan-20 


1.305 


BD69I 


2001-Feb-19 


1.374 


BD69J 


2001-Mar-16 


1.432 


BD69K 


2001-Apr-15 


1.501 


BD69L 


2001-May-19 


1.579 


BD69M 


2001-Jun-24 


1.662 


BD69N 


2001-Aug-16 


1.783 



Table 1. Observing dates and epochs. 




2.4514 2.4515 2.4516 2.4517 2.4518 2.4519 2.452 2.4521 2.4522 2.4523 
Julian Date . n e 



Figure 1. The light curve of R Cas provided by the AAVSO. The 
dashed vertical lines indicate the dates of our observations. 



3 RESULTS 

The optical light curve of R Cas provided by the AAVSC0, 
is shown in Fig. [1] The dashed lines indicate the VLBA 
observational epochs that span the time interval from 1999 
September 09 to 2001 August 16, corresponding to an optical 
phase interval of <f> — 0.158 to <f> = 1.783. Fig. [2] shows 
the flux density in Jansky of SiO maser emission toward R 
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Velocity (km/sec) 

Figure 2. The flux density detected by the KP antenna (solid 
line) and by the shortest baseline (LA - PT) (dashed line), for 
the two hands of polarization (red=right-hand, blue=left-hand). 



Epoch 


Radius (mas) 


error (mas) 


BD62A 


25.46 


0.67 


BD62B 


27.09 


0.05 


BD62C 


27.08 


1.77 


BD62D 


26.05 


1.22 


BD62E 


27.10 


1.59 


BD62F 


28.43 


1.98 


BD62G 


27.41 


2.11 


BD62H 


28.09 


0.3 


BD62I 


28.49 


0.53 


BD69A 


27.29 


1.64 


BD69B 


27.71 


0.37 


BD69C 


28.00 


1.5 


BD69D 


27.98 


0.16 


BD69E 


27.50 


0.5 


BD69F 


28.15 


1.71 


BD69G 


28.52 


2.5 


BD69H 


29.16 


0.26 


BD69I 


27.15 


0.17 


BD69J 


26.83 


0.63 


BD69K 


25.11 


3.11 


BD69L 


23.96 


1.22 


BD69M 


25.87 


3.83 


BD69N 


26.82 


1.68 



Table 2. Estimated radii of the SiO maser shell at each epoch. 



Cas in both RR and LL polarizations. In this figure, the 
solid line indicates the single-dish flux density for the Kitt 
Peak antenna, and the dashed line, the flux density for the 
shortest VLBA baseline PT-LA. It clear from the figure that 
the shortest baseline of the VLBA detected about 27 percent 
of the single dish SiO flux. Comparisons for all epochs are 
shown in Figs. Bl - B3. 

The gross morphology of the SiO maser emission toward 
R Cas in each individual epoch is shown in Fig.[3]and Fig. [4] 
Predominantly, the maser emission is confined to a narrow 
projected ring in phases cj> = 0.241 — 1.374. The first epoch, 
4> = 0.158, shows an additional, outer arc; a few outlying 
masers are also seen at (j> = 0.744 — 0.873. From <j> — 1.432 — 
1.783, the maser emission has an irregular shape, and the 
maser spots are very faint. The shell is dominated by an 
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Eastern arc during the first stellar period and by a Western 
arc during the second period. Fig. [5] shows the SiO maser 
features for all 23 epochs, superimposed on one plot. This 
shows that the SiO emission is predominantly situated in a 
ring that is poorly filled at the north and south. This ring 
has a radius of 1.6 — 2R+ depending on the obsering epoch. 

Fig. [6] shows the observed ring radii (tabulated in Ta- 
ble. [2} plotted as a function of stellar phase. The error bars 
on the points indicate the ring thickness, computed as the 
standard deviation in the observed ring radius. We fitted 
a parabolic function (solid line) to the data in Fig. [6] that 
shows the increase in the ring diameter over the p hase range 
4> ~ 0.1 to </> ~ 1.3 (|Diamond fc Kemballl [20031 ) . The pro- 
jected shell starts shrinking beyond the optical phase ~ 1.3. 



4 DISCUSSION 

4.1 Morphology 

The ring-like appearance of the SiO masers (Fig. [5} is 
thought to arise from a spherical shell, wherein strong radial 
acceleration or deceleration causes predominantly tangential 
amplification. The remnant of an outer, presumably older, 
shell is seen in the East in Fig. [3] (cj> = 0.158) and Fig. H At 
this first epoch and at <f> > 1.234, some faint emission is seen 
within 20 mas of the centre of expansion. This could be due 
to the emergence of a new shell and, indeed, the shell ra- 
dius shrinks at ^ 1.374 (discussed further in Section r4.3[) . 
Alternatively, the central masers could be radially beamed 
from parcels of gas temporarily moving at a fairly steady ve- 
locity, in between acceleration and deceleration. Most of the 
inner masers are red-shifted, suggesting infall, since they lie 
along the line of sight to the star, which would be optically 
thick to maser emission from the back of the shell. 

4.2 Maser shell asymmetry 

Figs. O and U show that the SiO emission is dominated by an 
Eastern arc during the first stellar cycle, whilst the West- 
ern side of the shell dominates for the remaining epochs. 
More than 3/4 of the emission comes from the East at 
4> = 0.158 — 0.744 and more than 3/4 comes from the West 
at (j) — 1-091 — 1.783 as it is seen in Fi g. [7] A similar domi- 
nance by one (Western) arc was seen bv lPhillips et al.l (|200ll ) 
in two epochs of observation of the 86-GHz SiO masers in 
1998. There are several possible explanations. 

Firstly, if conditions in one hemisphere of a spherical 
SiO maser shell are more favourable for masing, the change 
in E-W asymmetry from phase to phase could be due to 
bulk rotation, away from us in the E and towards us in the 
W, about an N-S axis. This would require an equatorial 
velocity of 10-20 km.s -1 which is about twice the observed 
line of sight velocities, but this could be due to turbulence 
and selective maser amplification in our direction. However, 
it is hard to conceive a physical mechanism for spherical 
solid-body rotation, although a flared disc might be possible. 

Secondly, the sta r may eject mass into ar cs, as has been 
seen in IRC+10216 (|Murakawa et. al.l 120051 ), which could 
mean that at some cycles, mass ejections were concentrated 
along the line of sight. In such cases we would see little or 
no SiO maser emission. 



Thirdly, ste l lar ac tivity may disrupt masing close to the 
star. iRudnitskiil l|2008l ) report that R Cas is known to display 
Ha flares at intervals of several years (also seen from other 
Miras). This is the most probable explanation although the 
second factor could play a part. 

4.3 Maser shell diameter 

Fig.[S]shows the variation of the ring diameter for 23 epochs 
as a function of stellar phase. It is evident that expansion 
is the dominant overall kinematic behaviour between opti- 
cal phases </> ~ 0.158 to <f> ~ 1.3. Thereafter, the diameter 
decreases until a phase of <j> ~ 1.58. This may be due to the 
disappearance of outer maser features and/or infall under 
the influence of gravity. At (j> > 1.58 the diameter increases 
again, suggesting the impact of a new stellar pulsation. Sim- 
ilar behaviour was seen in TX Cam (|Diamond fc Kemballl 
l2003ft . 

We found the gravitational acceleration (<?sio) was ~ 
3.27 x 10~ 7 km s -2 at the midpoint of SiO maser region 
(~ 26.56 mas) for a star mas s of 1.2 Mq. By adopti ng a 
distance to R Cas of 176 pc (|Vlemmings et all 12003! ) . we 
found that the required infall time for a parcel of gas across 
a maser zone of size ~12 mas is ~ 500 days. This estimate 
has an uncertainty of > 20 percent, mainly due to the un- 
certainty in the distance to R Cas, but is comparable to, 
or greater than, the stella r pulsation period of 431 days. 
iDiamond fc Kemballl (|2003h found a similar relationship for 
TX Cam. We also used our estimate of <?siO to calculate 
that, if a parcel of gas was stationary at <f> — 1.308, it would 
attain a velocity of 14 km s _1 at cj> = 1-579, which is almost 
equal to the velocity of 13.7 km s _1 required to reach the 
observed minimum radius. 

These results indicate that complex kinematics exist in 
the SiO maser region. The gas is subject to two main forces; 
outward pressure due to stellar pulsations and the inward 
force of gravity. Fig. [6] shows that expansion appears to be 
a slower process than infall over the same distance, con- 
sistent with a variable outward force but a more constant 
inward force. When the stellar pulsation pressure is weakest, 
at the outer edge of the SiO maser shell, the parcel of gas 
would start to fall back towards the star under its gravita- 
tional pull. If the parcel does not reach the inner edge of the 
shell before the next pulsation, it will then be pushed out 
futher. These results, like those for TX Cam, are a step to- 
wards quantative verification of pulsation models like lBowerJ 
|l98St ). 

4.4 Maser models 

We make some simple comparisons between the obser ved 
maser structure and the model of iGrav et. all (|2009l ); a 
more detailed analysis will be presented in a future paper. 
The model deals with tangentially beamed emission, but al- 
though emission is seen along the line of sight to the star 
at some epochs, this is too faint to influence comparison of 
the models with the average velocity profiles or radii. This 
model predicts that the diameter of the maser ring should 
be about twice the size of the stellar photosphere. The op- 
tical stellar diamete r of R Cas is 25.3 ± 3.3 mas at (j> = 0.93 
jWeigelt et-al-l feoOO). We found that the SiO maser mean an- 
gular diameter, averaged over 23 epochs, is ~ 54 mas which 
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is double the photospheric diameter within the uncertain- 
ties, showing that the model is consistent with our data. 

This model also predicts that the ring size increases 
from stellar phase 0.1 to 0.25 and then decreases. Fig. [6] 
shows that, during the first cycle, the radius continues in- 
creasing past phase 0.25. However, in the second cycle, the 
radius decreases sharply at the predicted phase.. The model 
predicts that the maser weakens between phases 0.25 and 
0.4. Fig. [S]shows that the flux density peaks at <f) ~ 0.2 and 
the then decreases over the same phase range which suggests 
good consistency between our observations and the model. 

Another prediction of the model is that the average ve- 
locity of the transition studied here moves to higher, red- 
shifted velocities, from phase 0.1 to 0.4. We tested this us- 
ing the single dish (autocorrelation) observations shown in 
Figs. Bl - B3, in order to ensure that the full spectra were 
measured. The red and blue, dashed lines in Fig.[9]show the 
velocity structure averaged over this range of phase for the 
first and second cycles, respectively. The black vertical line 
indicates the stellar velocity (with an uncertainty of 2 km 
s _1 ). The blue vertical line marks the flux weighted mean 
velocity for all epochs. The black, solid line is the model. 
This figure shows that the peak flux density is red-shifted 
in both the present R Cas data for the first cycle and in the 
model. The separation in the velocities of the peaks in the 
model and R Cas spectra is just 1.25 standard deviation in 
the R Cas data. The second cycle is not consistent with the 
model. 

Fig. [10] shows that the SiO maser variabilty in R Cas 
has approximately the same period as its optical light curve 
with a phase lag « 0.153 (66 d a ys) w hich is consistent with 
the predictions of iPardo et. al.l l|2004h . 

At this point, we compare the observational data of 
R Cas with a computational model of an AGB star. This 
model combines hydrodynamic sol utions with a mase r radia- 
tive transfer code as described in iGrav et. all (|2009h . apart 
from the exceptions d i scusse d below. 

The IGrav et. al.1 (|2009l ) model was extended to sam- 
ple 12 stellar phases between <f) = 0.4 to <f) = 1.5, so that 
the model pulsation passes through an optical maximum. 
We note that the m odel star is not spe cifically designed to 
represent R Cas. In IGrav et. al.l (|2009l ') . the 4 phase sam- 
ples were modelled as separate snapshots, unrelated to each 
other in the sense of positional evolution of the masing ob- 
jects. By contrast, the current model evolves the maser dis- 
tribution by taking the output condition s of one sample as 
the in put conditions for the next, as in iHumphrevs et. al.l 
(2002). Therefore the maser spots have random positions 
only for the first sample (at phase 0.4 in the first cycle). We 
computed the inner, outer and mean radius for each model 
epoch, defined in the same way as for the observational data. 
The mean radius of the maser shell is plotted for R Cas, in 
Fig. [TTJ It is evident that the model maser shell is substan- 
tially smaller than R Cas. However, a significant decrease in 
mean radius is apparent following optical maximum in both 
model and R Cas. This is the most statistically significant 
radius change in the observational data. This shrinkage is 
delayed in the R Cas data, with respect to the model, by 
about 0.2 stellar periods. On inspection of the model spec- 
tra during the radius drop, the peak of the model 43-GHz 
spectrum shifts to the red by ~7 km.s -1 . 
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Figure 3. Total intensity images [Epochs (1-9)]. The symbol size 
is proportional to the square root of the flux density of each com- 
ponent, the large circles represent the mean maser shell. 
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Figure 4. Same as Fig. [3] [Epochs (10-23)]. 
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Figure 5. SiO maser features observed at all 23 epochs. 




31 

30 

29 

_ 28 
£ 

9 27 

H 

I 25 
24 
23 
22 
21 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
Stellar phase 

Figure 6. The variation of the SiO maser ring radius with the 
stellar phase. 



This is consistent with the observed reddening of the 
R Cas emission during its contraction. We also inspected 
the model velocity at mean radius R, and noted that the 
velocity was closer to infall at the R during ring contraction 
than at optical maximum. 



5 CONCLUSIONS 

The SiO maser morphology in R Cas for the transition v 
I . J=l— > is described as a rarely completed ring structure, 
and the ring is dominated by an Eastern arc in the first 
stellar period and by a Western arc in the second stellar 
period. The result are compared with a model produced by 
iGrav et. al.l (|2009l ) which predicts most of observed features 
in our data. The infall time towards the star is greater than 
the pulsation period with a factor of ~ 1.2. The maser light 
curve behavior is similar to the optical light curve but with 
phase lag of 66 days. 
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Figure 7. The total East- (Red) and West- (Blue) flux density 
for 23 epochs. 
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Figure 8. The peak SiO maser flux density as a function of 
the stellar phase over the interval described by the model of 
IGrav et. al.l <200S|) . 
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Figure 9. Total intensity SiO maser velocity profiles measured by 
a single dish. The red and blue dashed curves show the emission 
averaged over the first and second cycles, respectively. The black 
solid curve represents the model. The black vertical line is the 
stellar velocity with an uncertainty (±2 km.s -1 ) and the blue 
vertical line is the flux- weighted mean velocity for all epochs. 
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Figure 10. Phase difference between the maser light curve (upper 
frame) and the optical light curve (lower frame). 



32 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 l.S 2.0 
Stellar phase 

Figure 11. A comparison between the behaviour of the radius 
in the model (blue) and our data (red). 
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APPENDIX A: ESTIMATES OF THE STELLAR 
VELOCITY 

Various estimates of the stellar velocity V+ appear in the 
literature; we tabulate these in Table |A"T1 and adopt the mean 
value, +24 km s _1 , standard deviation 2 km s~ . 
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Vlsr (km. s 1 ) 


Reference 


20.4 


Jewell et. al. (1991) 


24 


Colomer et. al. (20001 


26 


Vlemmings et al. (2003) 


25.4 


Boboltz et. al. (2004) 


24.9 


Matthews et. al. (2007) 



Table Al. Stellar velocity. 



APPENDIX B: PLOTS OF FLUX DENSITIES 
FROM A SINGLE DISH ANTENNA AND FROM 
THE SHORTEST BASELINE FOR DUAL-HAND 
POLARIZATION 

Figs. IB ll IB2I IB3l compare the flux density from a single dish 
Kitt Peak (KP) (solid line) and from the shortest baseline 
between Los Alamos and Pie Town (LA - PT) (dashed line) 
for the two hands of polarization (red = right hand, blue 
= left hand). The number on each plot indicates the stellar 
phase. The apparent offset between L and R is probably a 
residual error. 
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Figure Bl. The flux density for antenna KP (solid line) and shortest baseline (LA — PT) (dashed line) for the two hands of polarization 
(red=right-hand, blue=left-hand) for epochs 1-8. 
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Figure B2. Same as Fig. Bl for epochs 9 -16. 
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Figure B3. Same as Fig. Bl for epochs 17 — 23. 



